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ABSTRACT 

We investigate the detectability of the primordial CMB polarization B-mode power 
spectrum on large scales in the presence of instrumental noise and realistic foreground 
contamination. We have worked out a method to estimate the errors on component 
separation and to propagate them up to the power spectrum estimation. The perfor- 
mances of our method are illustrated by applying it to the instrumental specifications 
of the Planck satellite and to the proposed configuration for the next generation 
CMB polarization experiment COrE. We demonstrate that a proper component sep- 
aration step is required in order achieve the detection of B-modes on large scales and 
that the final sensitivity to B-modes of a given experiment is determined by a delicate 
balance between noise level and residual foregrounds, which depend on the set of fre- 
quencies exploited in the CMB reconstruction, on the signal-to-noise of each frequency 
map, and on our ability to correctly model the spectral behavior of the foreground 
components. We have produced a flexible software tool that allows the comparison 
of performances on B-mode detection of different instrumental specifications (choice 
of frequencies, noise level at each frequency, etc.) as well as of different proposed 
approaches to component separation. 



1 INTRODUCTION 



Primordial B-mode (curl component) polarization of the 
CMB provides a unique opportunity to detect the imprint 
of the primordial gravitational waves predicted by the in- 
flationary paradigm. Measuring the amplitude of these ten- 
sor perturbations would fix the energy scale of inflation and 
its potential and would provide a powerful constraint on a 
broad class of inflationary models. Moreover, the confirma- 
tion of inflation and the determination of the inflationary 
potential would have profound implications, by providing a 
direct observational link with the physics of the early uni- 
verse. 

The bulk of the statistical information on inflationary 
B-modes is concentrated in two features of the power spec- 
trum: the reionization bump at multipoles £ ~ 2 — 10 and 
the main bump at multipoles I ~ 20 — 100. Given that most 
of the signal lies on large angular scales (larger than ~ 1°), 
full-sky surveys, and thus satellite experiments, are required 
to limit cosmic variance. 

When designing a satellite experiment targeted to B- 
mode detection on large scales it is mandatory to consider, 
together with the signal-to-noise ratio, also the issues related 
to foreground contamination and component separation. In 
fact, as demonstrated by the WMAP data (Gold et al. 2009), 
the cleaning of CMB E-mode from polarized foregrounds 
is still a manageable problem. But polarized foregrounds 
are expected to dominate by a wide margin over the CMB 
B-mode because foreground E- and B-modes have similar 



amplitudes while the CMB B-mode is far weaker than the 
E-mode. 

The constraints in detecting primordial B-modes due 
to foreground contamination and residuals from foreground 
subtraction have been investigated by several authors. Many 
papers include a modeling of the residuals in the Fisher ma- 
trix (see. e.g., Tucci et al. 2005, Amarie et al. 2005, Verde 
et al. 2006) but generally bypass any realistic component 
separation approach and thus cannot properly deal with 
component separation errors. Other analyses focussed on 
a certain component separation method and evaluated its 
performances through simulations tailored for a specific ex- 
periment. For example, Betoule et al. (2009) and Efstathiou 
et al. (2009) propose different strategies (a spectral matching 
component separation method and an internal template fit- 
ting method, respectively) targeted to B-mode detection for 
the Planck experiment. More recently, Stivoli et al. (2010) 
assessed the capabilities of a parametric component separa- 
tion method to detect B-modes for suborbital experiments. 

The purpose of our work is different: while understand- 
ing the need of performing detailed simulations once the ex- 
perimental configuration is established, we want to provide 
a fast, agile and flexible toolQ to forecast the detectability of 
B-modes allowing us to easily change both the instrumental 
specifications and the CMB reconstruction strategy. Such a 
tool can be usually exploited to help optimizing the instru- 
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mental design and the data reduction pipeline for future 
CMB polarization experiments. 

The different components are reconstructed as a suit- 
able linear mixture of the data. Different choices for the 
linear mixture operator are available, and we carry out a 
comparison of their performances The forecast, targeted to 
large scales, on the B-mode detectability is performed at 
the power spectrum level and includes the computation of 
the final power spectra of noise and foreground residuals on 
the CMB map. We explicitly account for errors in the fre- 
quency scalings assumed for the foreground components. On 
the other hand, we do not account for errors due to actual 
power spectrum estimation (leakage effects) nor for instru- 
mental systematics. The former should be subdominant for 
the almost full-sky coverage considered here, especially when 
exploiting a power spectrum estimation method optimized 
for low multipoles (e.g. Gruppuso et al. 2009). Instrumental 
systematics obviously depend on the particular instrument 
and thus cannot be included in our general treatment. 

The outline of our paper is the following: in §[2] we de- 
scribe our models for the data and the component separation 
process. In §[3] we derive analytical expressions for the noise 
and residual foreground contributions to the errors in the 
BB power spectrum. In §|4]we apply our method to the spec- 
ifications of the Planck satellite and of a proposed CMB 
polarization experiment (COrE white paper in preparation). 
Our results are presented in §[5]and our main conclusions are 
summarized in §[U 



2 STATEMENT OF THE PROBLEM 
2.1 Data model 

The microwave sky contains, besides the CMB, several fore- 
ground components, both diffuse and compact. For our anal- 
ysis, which is focused on large and intermediate scales, we 
will consider only diffuse foregrounds. The main diffuse po- 
larized foregrounds are Galactic synchrotron and thermal 
dust (the free-free emission is unpolarized and the anoma- 
lous dust emission is also expected to be essentially unpo- 
larized) . 

The synchrotron component dominates at the lower fre- 
quencies. Its spectral behavior in antenna temperature can 
be modeled as a power law: 



TA,synch(^) OC V 



-0s 



(1) 



where the synchrotron spectral index p s can vary in the sky 
in the range 2.5 < f} a < 3.5 (Gold et al. 2009). The spectral 
behaviour of thermal dust emission, that takes over at high 
frequencies, follows approximately a grey-body law: 



?A,dust(f) OC 



exp(hv/kT dus 



(2) 



Both Pd and Td are spatially varying around /3d ~ 1-7 and 
Td ~ 18 K. The polarized CMB signal has a blackbody spec- 
trum: 



Taj 



t (u) OC 



(hu/kTcMB) 2 exp(hu/kT C MB) 
(exp(^/fcT CM B) - l) 2 



(3) 



with Tcmb = 2.73. 

The sky radiation, x, from the direction r at the fre- 



quency v is the superposition of signals coming from N c 
different physical processes Sji 



3=1 



(4) 



The instrument integrates the signal over frequency in each 
of its Nd channels, convolves it with a certain kernel, repre- 
senting the spatial response function, and adds noise. The 
measurement yielded by the generic channel centered at the 
frequency v, with beam axis in the direction f is: 

x v (f) — I B(r — r',v') t v (y')Sj{r', v')dr'dv + n v {f), (5) 

3=1 

where B(r — r',u') is the (azimuthally-symmetric) beam, 
tu(v') is the frequency response of the channel and n v (f) is 
the noise map. The data model in eq. ([5]) can be simplified 
under the following assumptions: 

(i) Each source signal, Sj(r,v), is a separable function of 
direction and frequency at least within limited sky patches; 

(ii) B{r,v) is independent of frequency within the pass- 
band of each channel; 

(iii) B(f , v) is the same for all the channels. 

The first two are reasonable approximations. To meet the 
third one we need to smooth the data to the lowest resolution 
we are dealing with. This is not a problem in the present 
context, since we are interested in large angular scales. 

Under these assumptions, the linear mixture data model 
applies. In vector form, for each direction of the sky (each 
pixel) we may write: 



x = Hs + n. 



(6) 



where x and n are A^-vectors containing data and in- 
strumental noise respectively, s is the iV c -vector containing 
sources, and H is the Nd x N c mixing matrix, containing the 
frequency scaling of the components. The spatial variability 
of the synchrotron and dust spectral indices implies that the 
mixing matrix H is in general different for different pixels. 



2.2 CMB recovery 

If the linear mixture data model holds, the components 
which are mixed in the channel maps x [eq. <j6j) ] can be 
reconstructed as: 



s = Wi. 



(7) 



where S is an estimate of the components s and W is a N c x 
Nd matrix called reconstruction matrix. We choose to rely on 
a linear estimator basically because it will allow us to easily 
deal with the component separation process in the forecast 
of B-mode detection, as we will see in the next section. In 
addition this approach is handy for Monte Carlo simulations 
needed to accurately control error sources. 

We adopt the so-called Generalized Least Square solu- 
tion (GLS): 



W = [H T N" 1 H]" 1 H T N _1 . 



(8) 



This requires the noise covariance N of the channel maps 
and an estimate H of the mixing matrix H that can be ob- 
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tained exploiting any of the successful component separa- 
tion methods discussed in the literature (see, e.g., Bonaldi 
et al. 2006, Eriksen et al. 2006, Stompor et al. 2009). We 
stress that this choice is not completely general, as other re- 
construction matrices could be used. The adaptation of our 
approach to other choices for W, still in the framework of 
the linear mixture data model, is however straightforward. 

The fact that our reconstruction matrix explicitly con- 
tains the estimated mixing matrix will allow us to easily 
introduce in our forecast also the effect of errors in the mix- 
ing matrix estimation (see §[3}. 

The noise covariance term, N, in eq. (JS| allow us to per- 
form the component separation keeping the noise level under 
control. Nevertheless, every component separation leads to a 
noise amplification by an amount ultimately depending on 
the conditioning of the mixing matrix H, as we will show 
in §31 Despite this drawback, in the low multipole regime 
considered we can not dispense from the component separa- 
tion step (e.g. by resorting to an aggressive masking of fore- 
ground contaminated regions) because the foreground con- 
tamination is more troublesome than the noise amplification 
and large areas are essential to restrain cosmic variance. 

On the contrary, the use of a suitable mask could be 
enough to keep foreground contamination under control 
on smaller scales, where the power spectra of diffuse fore- 
grounds decrease more rapidly than those of CMB and noise. 
In this case we could completely change strategy and sim- 
ply recover the CMB as a weighted mean of the channels 
with weights given by the inverse noise variance. This can 
be viewed again as a linear combination as in eq. (0, where 
the reconstruction matrix W has null elements for the fore- 
grounds and inverse noise variance weights for the CMB. 
This approach is complementary to the one previously de- 
scribed: it allows noise reduction but does not separate the 
components; only one component is reconstructed and con- 
sidered as CMB. 

In this work we will consider both approaches, which 
will be labeled as "component separation" (CS) and "mini- 
mum variance" (MV). The comparison of the results in §[5] 
will give us an idea of the multipole ranges where each ap- 
proach works better. 



3 FORECAST OF B-MODE DETECTION 

Our forecast is done at the power spectrum level and ac- 
counts for both noise errors and foreground residuals: 



AC I — ACe , noise + AClfon 



(9) 



In this section we formalize the estimate of the two contri- 
butions ACe. noise and ACe.ioicg to the error on the power 
spectrum, exploiting the model of the sky [eq. (HI}] and of 
the component separation process [eq. (0] described in the 
previous section. 

To compute the noise error AC^noisc we have to es- 
timate the noise bias on the CMB power spectrum, say 
JVcmb , i.e. the power spectrum of the noise in the recon- 
structed CMB map. This can be done with "noise-only" 
Monte Carlo simulations whereby several sets of realistic 
noise maps are generated and each of them is mixed with the 
matrix W. Then the power spectrum is computed for each 



simulation and, finally, all the noise power spectra are aver- 
aged. Here we make the simplifying assumptions of Gaussian 
white noise and spatially invariant foreground spectral prop- 
erties. This allows us to proceed analytically. In the frame- 
work of a realistic spatially-varying spectral model, we have 
verified that the analytic computation of the noise error with 
mean spectral dependencies over the sky is still enough ac- 
curate for the purpose of the forecast here considered. In 
fact, due to the intrinsic uncorrelation between noise and 
foregrounds, the mean level of the noise bias is predicted 
with good accuracy. 

The theoretical power spectrum of Gaussian white noise 
at the frequency v is: 



47rfsky 2 „2 



(10) 



where fsky is the sky fraction considered, Npix the number 
of pixels, <y% the noise variance at frequency v and B„(£) is 
the beam function applied to the channel maps to obtain a 
common resolution. As in our case we deal with polarization, 
a v is the noise of the detector at frequency v multiplied 
by \pl which, in the Gaussian white noise hipothesis here 
considered, gives the RMS of the noise in Q and U maps. 

Given the linearity of the CMB recovery process [eq. (0] 
and under the assumption of a spatially invariant recon- 
struction matrix, the noise bias is obtained by combining 
the channel noise spectra N '„ with the matrix W 2 : 



-ZVcMB = 22 W "'CMB N„, 



(11) 



where w„ cms are the elements of the matrix W 2 pertaining 
to the CMB component. 

The error on the CMB power spectrum is due to the 
fact that we do not know the actual noise realization, but 
only a mean noise bias. In other words, ACe , noise is the error 
due to the sampling variance of the noise bias ATcmb : 



AC 



£, noise 



2/(21 + 1) 
fsky nbin(£) 



Net, 



(12) 



where nbin is a function containing the number of multi- 
poles around any chosen I to be averaged according to a 
certain binning scheme. 

The error AC^forcg [eq. (©] due to the imperfect fore- 
ground subtraction can be computed following Stivoli et al. 
(2010). The map of residuals, s — s, for a linear mixture 
source reconstruction can be estimated as: 



8-8 = (WH - l)S, 



(13) 



where I is the identity matrix and § is a set of simulated 
components. ACf,f or eg is the power spectrum of the residu- 
als computed over a certain area of the sky. We note that 
in the computation of the foregorund residuals we do not 
proceed analytically up to the power spectrum level. This 
means that, in this case, we could easily exploit a realis- 
tic spatially-varying foreground spectral model. Given the 
lack of observational constraint at present, however, for the 
purpose of the paper we are relying again on the spatially- 
invariant model. 

According to eq. (|13[) the foreground residuals vanish 
for W = H _1 , which is the exact solution of the component 
separation problem in the absence of noise. For the recon- 
struction matrix of eq. J5}, the foreground contamination is 
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Table 1. Instrumental characteristics considered in the present study for the Planck and COrE experiments. 
The rms per pixel is referred to nside 1024 (size ~ 3.5arcmin) 



Planck specifications^ 


v (GHz) 


30 


44 


70 


100 


143 






217 








353 










FWHM (arcmin) 


3.3 


2d 


11 


9.5 


7.1 






5.0 








5.0 










29 m RMS AT(,uK RJ) 


66 


64 


60 


20 


10 






7.6 








7.5 










COrE specifications^ 


v (GHz) 




45 


75 


105 


135 


165 


195 


225 


255 


285 


315 


375 


435 


555 


675 


795 


FWHM (arcmin) 




2.3 


14 


10 


7.8 


6,1 


5.4 


4.7 


4.1 


3.7 


3.3 


2.8 


2.4 


1.9 


1.6 


1.3 


24 m RMS AT(fjK RJ)) 




1.4 


0.7 


0.6 


0.5 


0,1 


0,3 


0.2 


0.4 


0.5 


1.1 


3.6 


0.9 


1.2 


1.2 


1.3 



a LFI specifications as in Mandolesi et al. (2010), HFI specifications as in Lamarre et al. (2008) 
k ESA M3 call (Dec 2010) available at http://www.core-mission.org 



non null (although very low) even for H = H. It obviously 
increases with increasing error in the mixing matrix estima- 
tion. In the MV case the mismatch between W and H _1 is 
higher, and the foreground contamination increases. 

It is clear that the estimate of AC^f ore g through the 
previous relation is model-dependent since it relies on simu- 
lations of the data s which are hampered by our poor knowl- 
edge of polarized foregrounds. The situation will substan- 
tially improve in the near future as new polarization data, 
above all those from the Planck mission, will become avail- 
able. 



4 APPLICATION TO THE PLANCK AND 
CORE EXPERIMENTS 

The adopted experimental specifications for the Planck 
and COrE experiments are reported in Table 1. Planck 
specifications are derived from Mandolesi et al. (2010) for 
the Low Frequency Instrument (LFI) and from Lamarre et 
al. (2008) for the High Frequency Instrument (HFI) assum- 
ing a mission duration of 29 months. The specification for 
COrE are from the ESA M3 call (Dec 2010) available at 
http://www. core-mission, org. 

For each experiment we computed the error on the 
power spectrum using eqs. (f9"))- (|13l) for both the minimum 
variance (MV) and component separation (CS) case. 



gravitational lensing. We have added tensor modes with ten- 
sor to scalar ratio r=[0. 1,0. 03, 0.01, 0.001] and re-ionization 
optical depth r = 0.1. 



4.2 Minimum noise variance weighting 

As mentioned above, the MV approach assumes that, out- 
side a suitable mask, foregrounds are negligible so that we 
need to deal only with CMB and noise. Therefore the recon- 
struction matrix, W, contains only inverse noise variance 
weights for the CMB and null elements for the foreground 
components. For this assumption to be plausible we must 
restrict ourselves to a "CMB sensitive" set of channels and 
adopt a wide Galactic mask. To compute H and W we have 
thus considered only the channels in the frequency range 
70 ^ v < 200 GHz and computed the power spectrum of 
the foreground residual map [eq. (|13p ] only over the high 
Galactic latitude sky [|6| ^ 30°, so that fsky = 0.5, see 
eq. (H2J]. 

Given the different sensitivities of Planck and COrE 
we have used different bin sizes. Bin centers for I ^ 150 are 
tbin = [4, 9, 14, 20, 28, 39, 53, 80, 140] for Planck and l bin = 
[3, 6, 10, 19, 31, 49, 74, 106, 150] for COrE. 



4.3 Component separation weighting 

For the component separation approach we generated the 
estimated mixing matrix assuming the two standard po- 
larized foreground components (synchrotron and thermal 
dust) with spatially invariant spectra, but allowing for er- 
rors in the estimate of synchrotron and dust spectral indices. 
Three error regimes, labeled as "conservative", "intermedi- 
ate", and "goal" were considered (see Table [5}. 

The "conservative" errors on spectral indices are those 
estimated by Ricciardi et al. (2010) applying the CCA com- 
ponent separation method (Bonaldi et al. 2006) to a realistic 
simulation of Planck polarization data, including spatially- 
varying mixing matrix and realistic noise for 2 all-sky sur- 
veys. Thus they are slightly pessimistic for the real Planck 
mission, now extended to 4 surveys. Moreover the next gen- 
2 http://www.apc.univ-paris7.fr/APC_CS/Recherche/Adamis/PSM/pairyition experiments targeted to B-modes detection will un- 
en.php doubtedly improve the accuracy of the determination of fore- 



4.1 Sky model 

The polarization Q and U synchrotron and dust templates of 
our model were obtained at 100 GHz by running the Planck 
Sky Model (PSM).0 Extrapolations to lower and higher fre- 
quencies were made using the spectra of eq. ([TJ with /3 S = 3 
for synchrotron and of eq. ^ with fid — 1.7 and Td = 18 K 
for dust. Is has been necessary to slightly adjust the inten- 
sity of the synchrotron template to reproduce the WMAP 
7-yr K-band polarization map. The polarized CMB simu- 
lation is based on a standard ACDM model with WMAP 
7-yr cosmological parameters (Larson et al. 2010), including 



Forecast B-modes detection at large scales 5 



Table 2. Adopted errors in the estimation of the synchrotron and 
dust spectral indices for three different regimes 







Pd 


conservative 


0.05 


0.01 


intermediate 


0.01 


0.005 


goal 


0.005 


0.001 



ground spectral properties. That's why we also considered 
the "intermediate" and "goal" regimes. 

To propagate the errors on spectral indices to the final 
result, we generated a set of 10 estimated mixing matrices 
H, drawing the actual spectral indices from Gaussian distri- 
butions with a equal to the spectral index errors. For each H 
we computed the reconstruction matrix W, the residual map 
through eq. (|13p and its power spectrum. Since in this case 
we minimize foreground residuals the usable sky fraction is 
wider. We adopt fsky = 0.85, corresponding to a Galactic 
cut of ±10 deg around the equator. The values of AC^foreg 
shown in Fig. [T]are the mean for the 10 simulations. 

Another issue is the choice of the set of channels to be 
used for the source reconstruction. On one hand, using a 
wide frequency range brings in channels heavily foreground 
contaminated and, since the foreground spectra are not per- 
fectly known, this increases the foreground residuals in the 
final map. On the other hand, using a too limited set of 
channels in the linear combination leads to a higher noise 
level in the reconstructed CMB map, because the mixing 
matrix is less well-conditioned. 

The best trade-off clearly depends on the instrument, 
and in particular on its sensitivity. To exemplify this situ- 
ation we considered two sets of channels: a restricted one, 
limited to the frequency range 50 < v < 200 GHz, and a 
wider one (40 < v < 300 GHz) for both the considered 
experiments. 

We stress that this analysis on the choice of channels 
refers to the source reconstruction. For the estimation of the 
mixing matrix, i.e. of the foreground spectral properties, it 
is useful to exploit a frequency range as wide as possible. To 
this end the lowest and highest frequencies are particularly 
useful, as they map respectively polarized synchrotron and 
dust emissions with high signal-to-noise and low contam- 
ination from the other components. Even if not explicitly 
addressed in this paper, the mixing matrix estimation has 
a key role in the detection of B-modes, as will be shown by 
the comparison of the results on the spectral index errors 
for our different regimes. 



5 RESULTS 

In Fig. [1] we show the result of the forecast for Planck 
(upper panels) and COrE (lower panel) for the two choices 
of the frequency range (left: 50 < v < 200 GHz, right: 
40 < v < 300 GHz) used for the CMB map reconstruction. 
We can immediately note that the MV approach performs 
worse: the noise level is low, but the total error is strongly 
dominated by foregrounds even with a Galactic mask cov- 
ering 50% of the sky ^ 30°). To verify wether this re- 
sult depends critically on the mask, we also combined the 
|6| ^ 30° mask with the WMAP 7yrs polarization mask, fi- 



nally exploiting 46% of the sky. In this case, the foreground 
contamination lowers approximately by 30% at large scales. 
With a mask optimized for the frequency coverage here con- 
sidered we could do even better; nonetheless, the residual 
contamination would be still very high respect to the CMB. 
We conclude that on large and intermediate scales a proper 
component separation is necessary. 

In the Planck case, the total error for the CS approach 
is noise dominated even for a 85% sky coverage (|6| ^ 10°). 
Extending the frequency range from 50 < v < 200 GHz to 
the 40 < v < 300 GHz improves substantially the noise level. 
To be more quantitative, let's call A 2 the sum of the elements 
of the matrix W 2 pertaining to the CMB component. These 
elements weight the noise power spectra in the computation 
of the noise bias iVcMB- Thus A can be interpreted as the 
amplification of the noise rms. In our simulation, the smaller 
set of channels yields A ~ 3, while the larger one yields A ~ 
1.3. Our simulations indicate that, using the wider frequency 
range, Planck can detect B-modes for r = 0.1 and, with 
lower significance, for r = 0.03. The fact that the error is 
dominated by the noise term implies that the improvement 
in the mixing matrix estimation from the "conservative" to 
the "intermediate" case has a minor effect. 

Although our conclusions on Planck's capability for 
detecting B-modes are similar to those of Efstathiou & Grat- 
ton (2009), our approach differs substantially from theirs. 
These authors based their analysis on a minimum noise vari- 
ance combination of the 70, 100 and 143 GHz channels and 
assumed that foregrounds can be removed to high accuracy 
if the highest and lowest polarized Planck channels are 
used as foreground templates and the most contaminated 
areas (37% of the sky) are masked. However, as shown by 
our analysis, the foreground contamination corresponding 
to the minimum noise variance combination is very high 
even at high Galactic latitude. Therefore a sufficiently ac- 
curate foreground removal by template fitting is very chal- 
lenging. Conversely, a more standard component separation 
approach, as the one exploited for our forecast, keeps both 
foreground contamination and noise under control. 

The bottom panels of Fig. [T] show our expectations for 
the COrE mission. Due to the higher sensitivity of this ex- 
periment, the total error for the CS case is dominated by 
noise only for very low errors in the mixing matrix esti- 
mation ("goal"). For "conservative" errors in the spectral 
indices, the best results are obtained with the restricted set 
of channels. B-mode polarization can be detected down to 
r — 0.01 even with the "conservative" errors in the spec- 
tral indices. The minimum value of r that can be reached 
depends on our ability to estimate the mixing matrix. 



6 CONCLUSIONS 

We presented a method to forecast the detectability of 
CMB B-mode polarization on large and intermediate angu- 
lar scales given the experimental specifications of a full-sky 
mult i- frequency experiment. Our forecast accounts for both 
noise and residual foreground contamination after the CMB 
has been extracted from the data by means of a suitable 
linear combination of the frequency maps. The computa- 
tion of foreground residuals explicitly includes errors in the 
modeling of the data and in the estimation of the frequency 
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Figure 1. Forecast for B-modes detection by Planck (upper panels) and COrE (lower panel) for two choices of the frequency range 
used for the polarized CMB reconstruction: 50 < v < 200 GHz (left-hand panels) and 40 < v < 300 GHz (right-hand panels). Black solid 
lines: fiducial models for tensor to scalar ratios, r, of 0.1, 0.03 0.01 and 0.001. Dotted lines: cosmic variance limits (cyan for 50% sky 
coverage, grey for 85%). Coloured solid lines: total error for MV (red) and CS with conservative (blue), intermediate (green) and goal 
(magenta) errors in the mixing matrix estimations (see Table [2l . Those total errors are the sum of the noise (red triangles for MV, grey 
triangles for CS) and foreground (coloured squares) errors. 



scalings of the components. The forecast is quick and flexi- 
ble, and allows the selection of different sets of frequencies 
and of different methods for the CMB reconstruction. In 
this paper we considered the minimum noise variance (MV) 
and GLS component separation (CS) methods, but other 
methods can be easily implemented. 

We have investigated several issues: 

• How do the performances of minimum noise variance 
methods compare with methods exploiting component sep- 
aration? We have shown that, although component separa- 
tion leads to a noise amplification by an amount ultimately 
depending on the conditioning of the mixing matrix, in the 
low multipole regime considered in this paper we cannot dis- 
pense from the component separation step, e.g. by resorting 
to an aggressive masking of foreground contaminated re- 
gions, because the foreground contamination is more severe 
than that of noise and the mask cannot be too extended 
because large areas are essential to restrain cosmic variance. 

• Which is the optimal frequency range for the recon- 
struction of CMB polarization maps? Using a wide fre- 
quency range brings in channels heavily foreground contam- 
inated and, since the foreground spectra are not perfectly 
known, increases the foreground residuals in the final map. 
On the other hand, using a too limited set of channels in 



the linear combination leads to a higher noise level in the 
reconstructed CMB map, because the mixing matrix is less 
well-conditioned. We find that if the total error is dominated 
by noise, like in the case of Planck, it is convenient to use 
a rather broad frequency range (40 < v < 300 GHz) while 
in the case of more sensitive experiments, like the planned 
COrE mission, it is preferable to restrict ourselves to the 
"cleaner" frequency range 50 < v < 200 GHz. We stress, 
however, that this applies to the CMB reconstruction step. 
For the preliminary "mixing matrix estimation" step a much 
broader frequency coverage is essential. 

• How critical is a better understanding of polarized dif- 
fuse foregrounds for measurements of the power spectrum of 
primordial B-mode polarization? We find that in the case of 
Planck the main limitation comes from the detector noise. 
We estimate that with four all-sky surveys Planck can de- 
tect the re-ionization bump of the B-mode power spectrum 
for values of the tensor to scalar ratio r down to < 0.1, in 
agreement with the earlier conclusion by Efstathiou & Grat- 
ton (2009). On the other hand, to fully exploit the sensitivity 
of planned experiments specifically aimed at measuring the 
primordial B-mode anisotropies, it is essential to substan- 
tially improve the determination of spectral properties of 
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polarized foregrounds. Planck maps will allow an impor- 
tant step forward in this direction. 

As soon as more accurate multi-frequency maps of po- 
larized foreground emissions are available, our method will 
allow us to make more quantitative predictions that may 
help the design of next generation CMB polarization exper- 
iments and the optimization of the data analysis strategy. 
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